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Low-field magnetoresistance in p-type strained quantum wells is studied. It is shown that the 
Rashba mechanism leads to the cubic in quasimomentum spin-orbit splitting of the hole energy 
spectrum and the antilocalization behavior of low-field magnetoresistance is well described by the 
Hikami-Larkin-Nagaoka expression. 
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The combination of quantum coherence and spin ro- 
tation produces a number of interesting transport prop- 
erties. Numerous proposals for electronic devices that 
use spin-orbit coupling have appeared in last years, in- 
cluding gate-controlled sources and detectors of spin- 
polarized current Spin-orbit coupling results in the 
spin splitting of the energy spectrum when an inversion 
symmetry is lifted. The lack of inversion symmetry of 
the original crystal results in the splitting of the energy 
spectrum, which is linear and cubic in in-plane quasimo- 
mentum, k. This splitting is described by terms known 
as the Dresselhaus terms A In low-dimensional systems 
an additional mechanism of spin splitting is caused by 
the asymmetry of the confining potential (so called the 
Rashba terra?). In two-dimensional (2D) semiconduc- 
tor systems this asymmetry arises from asymmetry of 
the smooth electrostatic potential in the perpendicular to 
the 2D plane direction, from Schottky barrier potential, 
from asymmetry in doping layers disposition, and com- 
position gradient along the growth direction. It is very 
important that this asymmetry can be controlled by gate 
voltage. For electron 2D states, the Rashba term is linear 
in k. For hole 2D systems, the situation becomes more 
complicated because of four-fold degeneracy of the top- 
most valence band Tg of the parent material. Theoretical 
considerations of this problem and experimental studies 
show that the splitting is cubic in k in this case 

The measurements of interference induced low-filed 
magnetoresistance are the powerful tool for studies of the 
spin-splitting, spin- and phase- relaxation mechanisms. 
At present, there are numerous studies of n-type 2D 
systems2ii2ii±ii2iiiii2ii£iiS, whereas the more complicated 
p-type systems are studied noticeable ! e s9iLi§ii2i22i2i 
(for more references see review article by Zawadzki and 
Pfeffer— ). As for the strained quantum well, the antilo- 
calization and spin relaxation in 2D hole gas are practi- 
cally not investigated in these systems. 

In this paper, we present results of experimental study 
of the low-field magnetoresistance caused by spin relax- 



ation in p-type strained GaAs/In^Gai-ajAs/GaAs quan- 
tum well structures. It has been found that the mag- 
netoresistance shape is well described by the Hikami- 
Larkin-Nagaoka (HLN) expression 2 ^ that means that the 
leading term in the spectrum splitting is cubic in quasi- 
momentum. We show that in contrast to n-type systems, 
where such a finding implies that the Dresselhaus spin 
splitting mechanism is the main, the Rashba mechanism 
is responsible for the spin splitting of the hole energy 
spectrum in strained quantum wells under investigation. 

The GaAs/In^Gai-^As/GaAs heterostructures were 
grown by metal-organic vapor phase epitaxy on semi- 
insulator GaAs substrate. The quantum well was biax- 
ially compressed due to the lattice mismatch between 
In^Gai-^As and GaAs. Two types of heterostructures 
were studied. The structures of the first type, 3855, 
3856, and 3857, consist of a 250 nm-thick undoped GaAs 
buffer layer, carbon 5-layer, a 7 nm spacer of undoped 
GaAs, a 8 nm In0.2Ga0.sAs well, a 7 nm spacer of un- 
doped GaAs, a carbon (5-layer and 200 nm cap layer of 
undoped GaAs (see Fig.^l. The structure of the second 
type, 3951, was analogous, the only difference was the 
wider spacer, 15 nm, and as sequence the higher mobil- 
ity. The structures within the first group differ by car- 
bon density in <5-layers. The parameters of the structures 
are presented in Tabled The samples were mesa etched 
into standard Hall bars and then an Al gate electrode 
was deposited by thermal evaporation onto the cap layer 
through a mask. Varying the gate voltage V g from —IV 
to +3 V we changed the hole density in the quantum 
well from 3 x 10 11 cm~ 2 to 1 x 10 12 cm -2 . The analy- 
sis of the temperature dependence of the Shubnikov-de 
Haas oscillations showed that the hole effective mass was 
equal to (0.160 ± 0.005) m and did not depend on the 
hole density. 

Figure |3 shows the low-field magnetoconductivity, 
Aa(B) = p-^(B) - p^(0), measured at T = 0.44 K 
for structure 3857 as a function of a normalized magnetic 
field b = B/B tr , where B tr = H/ (2el 2 ) with I as the mean 
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FIG. 1: The cross section of structure 3857 and its energy 
diagram. 
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free path, is presented in Fig. [2J The antilocalization 
maximum at B = in the conductivity-versus-magnetic 
field curves decreases with lowering hole density and dis- 
appears at V g = 2.75 V, when p ~ 3.8 x 10 11 cm -2 . In 
the structures 3951 with the higher hole mobility the 
disappearance happens at p ~ 3 x 10 11 cm~ 2 . 

Theoretically, the low-fi eld ano malous magnetoresis- 
tance was studied in Refs. llll23l25l It was shown that 
when the spin splitting is cubic in k, the magnetoconduc- 
tivity curve should be described by the Hikami-Larkin- 
Nagaoka (HLN) expression 
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FIG. 2: The magnetoconductivity plotted against the re- 
duced magnetic field for different gate voltages, structure 
3857, T = 0.44 K. 



where W (9) stands for the probability of the scattering 
by an angle 9. 

For the Dyakonov-Perel spin-relaxation mechanism, 24 
which is dominant at low temperatures, the value of t s 
is determined by the spin-orbit splitting of the energy 
spectrum, MI3 oc k 3 , as follows 



— = 2f2 2 T 3 

T, 



(3) 



where T3 is defined by Eq. J2J. 

Taking into account both the cubic and linear terms 
leads to more complicated expression which was obtained 
(l)in Ref.l25l The following two parameters describing the 
spin relaxation arise in this case 



Here, Go = e 2 /(2n 2 h), and r s are the phase and spin 
relaxation times, respectively, ip(x) is the digamma func- 
tion, and r ra , n = 1, is the transport relaxation time 



1 



= / W(0)(1 - cosn 



(2) 



and 



= 2ft 2 n 

T 



= 2{nj n + ft 2 T 3 ), 



(4) 



(•5) 



TABLE I: The parameters of structures investigated 



Structure 


Ni (cm" 2 )" 


iV 2 (cm" 2 ; 


a p (cm -2 ) n 


(cm 2 /Vs) 


3855 


4 x 10 11 


3 x 10 11 


4.7 x 10 11 


4800 


3856 


8 x 10 11 


6 x 10 11 


7.5 x 10 11 


5700 


3857 


1.2 x 10 12 


8 x 10 11 


9.5 x 10 11 


8000 


3951 


1.2 x 10 12 


8 x 10 11 


5.4 x 10 11 


13100 



a N\ and N2 are the carbon density in outer and inner 5-layers, 
respectively. 



where Ml\ is the linear in k, Ml\ cx k, spin-orbit splitting. 

Comparison of the experimental data with theoreti- 
cal expressions for two limiting cases, when only the 
cubic or linear term is taken into account, is shown in 
Fig. 13 To span the characteristic minimuma in Act- 
versus-B curves, the fitting interval has been chosen as 
— 0.3 B tr < B < 0.3 B tr . Strictly speaking, the bound- 
aries of this interval do not satisfy the diffusion approx- 
imation B <C B tr in which framework the formulae for 
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FIG. 3: The magnetoconductivity as a function of normalized 
magnetic field for the two gate voltages: V g = 1.5 V (p = 
8x 10 11 cm -2 , ti = 5.4 x 10" 13 s) and V g = 2.25 V (p = 5.8 x 
10 11 cm" 2 , n = 3.9 x 1(T 13 s); structure 3857, T = 0.44 K. 
Symbols are the experimental data. The dotted lines are the 
best fit by the expression from Ref. |^j when only the linear 
in k term is taken into account. Dashed lines are the best fit 
by the HLN expression, Eq. Solid lines are the results 

of the simulation procedure (see Appendix EJ which is valid 
beyond the diffusion approximation. The fit has been done 
within the magnetic field range —0.3B tr < B < 0.3B tr - The 
fitting parameters are given in Table ITT1 



magnetoconductanc o 23 i 25 has been derived. Neverthe- 
less, one can see that the taking into account only the 
linear term does not allow us to describe satisfactorily 
the magnetoconductivity shape within the fitting inter- 
val while the HLN expression gives a good agreement. 
Beyond the diffusion regime, the HLN theory was gener- 
alized by Zduniak et alm& However, the final expressions 
are very complicated and inconvenient to use in the fit- 
ting procedure. Because of this, we used the simulation 
approach described in our paper, Ref. [2(J To take into 



TABLE II: The parameters of the best fit for the data pre- 
sented in Fig. |3 as obtained taking into account only the lin- 
ear in k term*— only the cubic in k term in the diffusion 
approximation;— and the cubic in k term beyond the diffu- 
sion approximation (see Appendix [XJ. 



Vg (V) 


Theory 


Tl/T0 


Tl/T a 




Ref. 25. 


0.020 


0.178 


1.5 


Ref. 23 


0.016 


0.051 




Appendix El 


0.014 


0.040 




Ref. 25. 


0.034 


0.142 


2.25 


Ref. 23 


0.017 


0.032 




Appendix lAl 


0.013 


0.025 
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FIG. 4: The temperature dependence of the phase and spin 
relaxation time as obtained from the fit of the experimental 
data by Eq. Q for structure 3857 at V g = 1.5 V. Solid line is 
the T _1 -law. 



account the spin-relaxation processes, Eq. (20) from this 
paper was modified as described in Appendix 1X1 (details 
will be published elsewhere). As Fig. illustrates, with 
the use of Eq. QAlj l we obtain a nice coincidence of calcu- 
lated and measured curves over the whole magnetic field 
range. Although the theory beyond the diffusion approxi- 
mation describes the magnetoresistance curve better, the 
fitting parameters ti/t^ and t\/t s are found to be close 
to those obtained with the help of Eq. (JIJ (see Table [n}. 
Therefore, it seems natural to analyze the experimental 
results with the use of the more simple HLN expression. 

Further indication, that just the cubic in k splitting 
is responsible for the spin relaxation, is reasonable be- 
havior of the fitting parameters obtained from Eq. 
with the temperature change. As seen from Fig.0] the 
parameter exhibits the behavior close to T _1 -law that 
corresponds to the phase relaxation caused by inelastic- 
ity of electron-electron interaction^ The parameter t s 
is temperature independent as should be for degenerated 
electron gas. Such analysis has been carried out for all 
the structures investigated and the results are collected 
in Fig.0 as a plot of the spin relaxation time r s against 
the hole density controlled by the gate voltage. 

For the first sight the fact that the magneto conduc- 
tance curves are well described by the HLN expression 
means that the Dresselhaus cubic term gives the main 
contribution to the spin splitting in the structure inves- 
tigated. Whether or not it is so, one can understand an- 
alyzing the hole density dependence of spin-orbital split- 
ting, Mls(p). For the Dresselhaus mechanism, the split- 
ting should be proportional to p 3 ^ 2 because O3 = 7fc 3 /4, 
where 7 is constant depending only on the band param- 
eters of the parent material (see Appendix A in Ref. ^] 



4 




5.0x10 



1.0x10 



P (cm" ) 



FIG. 5: The spin relaxation time as a function of hole density 
controlled by the gate voltage for all structures investigated, 
T = 0.44 K. 
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FIG. 6: The hole density dependence of to (open symbols) 
and n (solid symbols) . Solid lines are provided as a guide for 
the eye. 



for details). Experimentally, the value of spin-orbit split- 
ting M^3 can be found from Eq. using r s obtained 
above and T3. How the quantity T3 has been obtained is 
considered below. 

As seen from Eq. J5J) the relaxation time T3 is de- 
termined by the scattering anisotropy via the function 
W(9). Just the same function determines the rela- 
tionship between the quantum and transport relaxation 
times, To and t%, respectively. Therefore, we estimate T3 
using the experimental value for To 

t 3 - t y W{6)d6 1 J W{9){\ -cos30)d0, (6) 

and conceiving the physically reasonable angle depen- 
dence for W{9) so that the ratio 

K i = JW(6)(1- cos 6)d6/ Jw(9)d9, (7) 

to be equal to the experimental quantity tq/t\. The value 
of To has been obtained from the analysis of the magnetic 
field dependence of the amplitude of the Shubnikov-de 
Haas oscillations, while T\ has been found from the mo- 
bility value, ti = fj,m/e. The experimental hole density 
dependences of To and t\, presented in Fig-IHI show that 
the scattering is really anisotropic in all the structures 
and the to to t\ ratio lies in the interval from 0.2 to 0.5. 
This seemingly points to the fact that the scattering is 
mainly determined by ionized impurities and W(9) can 
be chosen in the form obtained, e.g., in Ref.l29l How- 
ever, our estimation shows that the electron mobility in 
this case should be one-two order of magnitude higher 
than that observed experimentally. We suppose that the 



roughness of the quantum well interfaces restricts the mo- 
bility in our structures. This mechanism is theoretically 
studied in Ref. |3(], where the explicit form for W(9) is 
derived and it is shown that the scattering anisotropy 
strongly depends on the parameter A characterizing the 
fluctuations of the quantum well width due to interfaces 
roughness. Using the form for W{9) from this paper we 
have chosen such a value of parameter A which satisfies 
the equality between the experimental value of To /ti and 
the calculated from Eq. Q value of Kqi- Then, with this 
A-value we have calculated T3 to to ratio. Doing so we 
have found that T3 = (0.7 . . . 0.8)to when i^oi lies within 
actual for our case range, K 01 = 0.2 .. . 0.5. Note, the 
close results, t 3 ~ t , are obtained if one uses W{9) cor- 
responding to scattering by remote ionized impurities^ 

Finally, we arrive the key figure of the paper, Fig.0 
where the value of spin splitting = h/y/2r^T s is plot- 
ted as a function of the hole density. One can see that 
(i) we do not observe the characteristic for the cubic in 
k spin-orbit splitting p 3 / 2 -dependence and (ii) different 
structures demonstrate significantly different value of the 
splitting for a given hole density. Both these facts unam- 
biguously show that no the Dresselhaus mechanism is 
responsible for the cubic in k spin-orbit splitting of the 
hole spectrum. 

Let us now discuss specific features of the Rashba ef- 
fect for holes in strained quantum well heterostructures. 
In general, this effect in hole 2D systems is considered in 
Refs. l6l7U8l l9l. Below we write out only the main expres- 
sions which help us to describe the experimental results 
presented quantitatively. We restrict our consideration 
by the case when only three hole bands are taken into 
account. They are the heavy- and light-hole Tg bands 



5 



and split off by spin-orbit interaction IV hole band. In 
this case, the energy spectrum is described by the 6x6 
Luttinger-Kohn Hamiltonian 31 which includes the terms 
responsible for the strain^ As shown in Ref.|U the 6x6 
Hamiltonian can be decoupled into two independent 3x3 
matrices of the form 



where ip and E are solutions of the Schrodinger equation 



C±iB FTi\/t/2B 

V2TiB/V2 F±iy/3/2B D J 



H = 

where 

A ± = 
B = 
C = 
D = 
F = 



Here, 7j stand for h 2 j^/(2mo), where 7^ are the Lut- 
tinger parameters, k z is the wave vector along the [001] 
growth direction, k 2 = k 2 + k 2 , 9 is the angle between 
the in-plane wave vector and the [100] direction, V(z) is 
the macroscopic electric potential in the heterostructure, 
Er s and E Tl are the energies of edges of corresponding 
bands, and 



(8) 



-(71 T 2 72 )fc' - (71 ± I2)k 2 + E T8 (z) + V(z) ± S 
2V3"/ 3 kk z 

V3fc 2 ( 7 | cos 2 26» + 7! sin 2 2(9) 1/2 
-^(k 2 z +k 2 ) + Er 7 {z) + V(z) 
2 72 (V2fc 2 -k 2 /V2). 



S 



a + l\ Aa 
a - 1 



(9) 



is the splitting of the Tg band due to strain caused by 
the lattice mismatch between GaAs and In^Gai-^As. In 
Eq. ©, b stands for the axial deformation potential of 
the valence band, a is the Poisson's ratio, Aa is the lat- 
tice mismatch between materials of the quantum well and 
barrier, and a is the lattice constant of the quantum well 
material. Let us estimate characteristic energies for the 
case of GaAs/Ino^Gao.sAs heterostructure. The value of 
Aa/a is about 1.4 %, a is approximately equal to 1/3, b 
is about —1.7 eV so that the value of strain induced split- 
ting 2\S\ is approximately equal to 90 meV. This value is 
five-ten times greater than the Fermi energy in our case. 
We find the Rashba splitting of the hole energy spectrum 
using the ratio Ep/(2S) as a small parameter. The band 
parameters 7^ and A = Er 7 — Ep s are supposed indepen- 
dent of z-coordinate. Then, in isotropic approximation, 
72 = 73 =7, the energy spectrum of the upper split off 
band for our case can be written as follows 



e ± ~ e ± m* 



(10) 



with 
M2 3 



6 7 2 fc 3 J dz\^\ 2 ^ 



1 



1 



E-E T& {z)~A-V(z) 



E-E Vs (z)-S-V(z) 

(11) 



A + ip = Eip. 



(12) 



It is clearly seen from Eq. (|llfl , that the Rashba splitting 
for all 2D subbands formed from the upper hole band is 
cubic in k in contrast to the electron energy spectrum 
where it is linear in k. 

Now we are in position to compare the experimen- 
tal ^.f^-versus-p dependences with those calculated from 
Eqs. (fTTfl . Q . To find the electric potential V(z), the 
Schrodinger equation has been self-consistently solved 
with the Poisson equation. We have used the param- 
eters jf for Ino.2Gao.8As and the value of band offset 
SE V — i?r s (GaAs) — £?r 8 (Ino.2 Gao.sAs), which are ob- 
tained by the linear interpolation from the values of GaAs 
and InAs: 7^ = -9.5, 7 L = -3.5, SE V = 75 meV, 
A = 0.35 eV (the signs of the Luttinger parameters corre- 
spond to the increasing of energy into the valence band) . 
As an example, the energy profile and the wave function 
for structure 3857, V g — 0, is shown in Fig. ^ To de- 
scribe the experimental f^-versus-p dependence for each 
structure, the parameter S has been used as a fitting 
one. One can see from Fig. that we are able to de- 
scribe well the experimental results obtained for different 
samples in our model only slightly varying the S'-value 
from the one to other structure. The different value of 
strain induced splitting for different samples seems to be 
natural. It can result, for instance, from deviation of 
In-content from its nominal value. As for the value of 
the strain induced splitting, 2\S\ ~ 75 — 90 meV, it cor- 
responds to the lattice mismatch and In-content laying 
within the intervals (1.2 - 1.4)% and (17 - 20)%, respec- 
tively. Let us direct attention to the interesting detail. 
The fif^-versus-p dependence exhibits the correspond- 
ing to Eq. (|llfl behavior, /1SI3 cx p 3 ^ 2 , only at low hole 
density, p < 2 x 10 11 cm~ 2 . At higher hole density this 
dependence has a maximum and sign change (not shown 
in figure). This feature is caused by the fact that the 
hole density is varied by means of variation of the gate 
voltage. Applying the gate voltage we change not only 
the value of the Fermi quasimomentum but the energy 
profile of the quantum well as well. In this case the in- 
tegral in Eq. (|llf> is not constant any more and gives 
additional p-dependence in Q 3 . Vanishing of spin-orbit 
splitting at some hole density means that the quantum 
well in this point becomes effectively symmetric. We real- 
ize that the approximations of large strain induced split- 
ting and z-independence of 7$ parameters made above 
are crude enough. Moreover, the well boundaries can be 
smooth and different, and the In-content can vary across 
the quantum well. These factors being taken into account 
could in principle change the value of S obtained from the 
fit. However, this should not change our interpretation 
in the large. 

In summary, we have shown that the Rashba mecha- 
nism results in the cubic in k spin-orbit splitting of the 
hole energy spectrum in strained heterostructures. The 
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FIG. 7: The hole density dependence of the spin-orbital split- 
ting for different samples. Symbols are the experimental data 
obtained as hQ^ — h/\/2T3T s , solid lines are calculation results 
(see text), dashed lines show p 3 ^ 2 -law for structures 3856 and 
3951. In brackets, the values of the fitting parameter S for 
each structure are shown. 



magnetoresistance curve in this case is well described by 
the HLN-expression, that allows us to find the spin split- 
ting as a function of the hole density. We have found that 
these dependence is nonmonotonic at relatively high hole 
density due to sensitivity of quantum well profile to the 
gate voltage. 
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APPENDIX A: NUMERICAL SIMULATION OF 
ANTILOCALIZATION 



The weak localization phenomenon for spin-less par- 
ticle using the numerical simulation of a particle mo- 
tion over the plane with randomly distributed scatter- 
ers is studied both within and beyond diffusion regime in 
Ref. |2(J It has been shown that obtaining from the simu- 
lation procedure the parameters of closed paths, one can 
calculate the quantum interference correction to the con- 
ductivity and its magnetic field dependence (see Eq. (20) 
in the paper cited). Taking into account the spin re- 
laxation processes leads to the following expression for 
the interference quantum correction (in more detail this 
generalization will be considered elsewhere)^ 



Sa(b) = - 27rG ° cos(bSi) exp (-fcy) 



-~ + exp (-1^) + i exp {-2la s ) 



,(A1) 



where I s is the total number of paths, d is the diameter 
of the area from which the particle starts to walk and in 
which it returns, k and Si are the length and algebraic 
aria of j-th closed path, 7 and 7 S are the phenomenolog- 
ical parameters describing the phase and spin relaxation 
and corresponding in real systems to ratios ti/t^ and 
ti/t s , respectively, the lengths and areas in this expres- 
sion are measured in units of mean free path and squared 
mean free path, respectively, and summation runs over 
all closed paths. In order to treat the experimental re- 
sults presented in this paper, we have firstly collected 
the parameters of closed paths U and ^ simulating the 
motion of particle as described in Ref. 1261 and, then, we 
have used Eq. I|A1[) to fit the experimental curves with 7 
and 7s as the fitting parameters. 
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